Scaling properties of discontinuous maps 
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We study the scaling properties of discontinuous maps by analyzing the average value of the 
squared action variable I 2 . We focus our study on two dynamical regimes separated by the critical 
value K c of the control parameter K: the slow diffusion (K < K c ) and the quasilinear diffusion 
(K > K c ) regimes. We found that the scaling of I 2 for discontinuous maps when K <C K c and 
K ^> K c obeys the same scaling laws, in the appropriate limits, than Chirikov's standard map in 
the regimes of weak and strong nonlinearity, respectively. However, due to absence of KAM tori, 
we observed in both regimes that I 2 oc nK 13 for ti>l (being n the n-th iteration of the map) with 
/3 » 5/2 when K < K c and (3 w 2 for K > K c . 

PACS numbers: 05.45.-a, 05.45.Pq 



I. INTRODUCTION AND MODEL 

Chirikov's standard map (CSM), introduced in Rcf. [1], 
is an area preserving two-dimensional (2D) map for ac- 
tion and angle variables (I, 9): 



I n+1 = I n + Kf(6 n ) , 

9 n +i = n + I n+1 , mod - 2n , 



(1) 



where f(9 n ) = sm(9 n ) (due to this choice of /(#), CSM is 
identified as a continuous map). CSM describes the situ- 
ation when nonlinear resonances are equidistant in phase 
space which corresponds to a local description of dynam- 
ical chaos [2]. Due to this property various dynamical 
systems and maps can be locally reduced to map (1). 
Thus, CSM describes the universal and generic behav- 
ior of nearly-integrable Hamiltonian systems with two 
degrees of freedom having a divided phase space com- 
posed of stochastic motion bounded by invariant tori 
(also known as KAM scenario) [2]. 

CSM develops two dynamical regimes separated by the 
critical parameter K c [1—7]. When K < K c , regime of 
weak nonlinearity, the motion is mainly regular with re- 
gions of stocasticity and I is bounded by KAM surfaces. 
See for example Fig. 1(a) where we present the Poincare 
surface of section for CSM with K = 0.01. Here, the 
value of K is so small that the Poincare surface of section 
is equivalent to the phase portrait of a one-dimensional 
pendulum. At K = K c , the last KAM curve is de- 
stroyed and the transition to global stocasticity takes 
place. Then, for K > K c , regime of strong nonlinear- 
ity, / becomes unbounded and increases diffusively. See 
for example the Poincare map of Fig. 1(b) where a single 
trajectory has been iterated 3 x 10 4 times. 

Even though CSM describes the universal behavior of 
area-preserving continuous maps, another class of Hamil- 
tonian dynamical systems is represented by the discon- 
tinuous map [8] : 



where f(9 n ) — sin(0 n )sgn[cos(0 n )]. Examples of physical 
systems described by discontinuous maps are 2D billiard 
models like the stadium billiard [9, 10] and polygonal 
billiards [11, 12]. The origin of the discontinuity in map 
(2) are the sudden translations of the action under the 
system dynamics. 

As well as CSM, map (2) is known to have two differ- 
ent dynamical regimes, however both diffusive, delimited 
by the critical value K c = 1/T [8] . The regimes K < K c 
and K > K c are known as slow diffusion and quasilin- 
ear diffusion regimes, respectively. As an example of the 
dynamics of map (2), in Fig. 2 we show typical Poincare 
surface of sections in both regimes (for comparison pur- 
poses we have used the same values of K as in Fig. 1 
for CSM). On the one hand, as can be observed by con- 
trasting Figs. 1(a) and 2(a), the main difference between 
CSM and map (2) is that for K < K c the later does not 
show regular behavior. In fact, due to the discontinuities 
of f{9) in map (2), KAM theorem is not satisfied and 
map (2) does not develop the KAM scenario. Since for 
any K ^ the dynamics of map (2) is diffusive, a single 
trajectory can explore the entire phase space. However, 
in the slow diffusion regime the dynamics is far from be- 
ing stochastic due to the sticking of trajectories along 
cantori (fragments of KAM invariant tori), see Fig. 2(a). 
On the other hand, for K > K c map (2) shows diffusion 
similar to that of CSM, compare Figs. 1(b) and 2(b). We 
want to add that independently of the value of K ^ 0, 
map (2) has five period-one fixed points at / = and 
9 = [0,7r/2,7r,37r/2,27r]. 

In particular, in Ref. [13] a scaling analysis of CSM 
was performed by studying the average value of I 2 as a 
function of K and the n-th iteration of the map. There, 
the following scaling law was reported: 



P 



(3) 



In+1 
9 n +l 



In + Kf(9 n ) , 

9 n + TI n+1 , mod - 2tt 



(2) 



where a « 2 for K <§; K c and small n while a « 1 for 
K ^> K c and large n, with j3 ~ 2 in both cases. The 
scaling (3) has also been validated for several dynam- 
ical systems represented by the standard map such as 
the Fermi- Ulam model [14-18], time-dependent poten- 
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FIG. 1: Poincare surface of section for CSM, Eq. (1), with (a) 
K = 0.01 and (b) K = 10. In (a) 20 initial conditions with 
9o = 7r and /o = [—0.5, 0.5) were iterated 10 3 times. In (b) a 
single initial condition with 6o = 3 and Io = 0.01 was iterated 
3 x 10 4 times. 



FIG. 2: Poincare surface of section for the discontinuous map 
of Eq. (2) with (a) K = 0.01 and (b) K = 10. T = 1. A 
single initial condition with 9o = 3 and Io = 0.01 was iterated 
3 x 10 4 times. 



tial wells [19], and waveguide billiards [18, 20]; among 
others [21, 22]. 

Although map (2) has the same structure than CSM, 
a systematic investigation of the scaling properties of 
I 2 for discontinuous maps is not available in the liter- 
ature. Thus, in this paper we undertake this task. For 
this purpose, here we study the properties of the map 
of Eq. (2) by analyzing the scaling of the average value 
of the squared action variable I 2 as a function of n, K , 
and Jo- We choose the scaling approach to I 2 reported 
in Ref. [13] because of the similarity of maps (1) and (2). 
Moreover, since map (2) shows diffusion in both dynami- 
cal regimes (K < K c and K > K c ), we expect the scaling 
(3) to be valid for discontinuous maps when diffusion is 
present with scaling exponents a and (3 to be determined. 



II. RESULTS 

We compute I 2 for map (2) following two steps [13]: 
First we calculate the average squared action over the 
orbit associated with the initial condition j as 

1 " 

( l2 n,j) = ^E4- ' 
i=0 

where i refers to the i-th iteration of the map. Then, the 
average value of I 2 is defined as the average over M inde- 
pendent realizations of the map (by randomly choosing 



values of 0o)' 

1 M 

I 2 {n,K,I ) = -Y J {ll, j ) ■ (4) 
In the following, without lost of generality, we set T = 1. 

A. Slow diffusion regime 

In Fig. 3(a) we present I 2 as a function of n in the slow 
diffusion regime (K cl) for several combinations of K 
and Iq- In fact, I 2 is always an increasing function of n, 
however its growth is marginal in some iteration intervals 
producing plateaus in the curves I 2 vs n. 

For / <c K, see full symbols in Fig. 3(a), I 2 grows 
up to a crossover iteration number n^. When n^) < 
n < n[ 2 ^ the trajectories wander around the period-one 
fixed points of the map making the growth of I 2 negli- 
gible; that is, I 2 becomes almost constant. We call this 
constant J 2 at . Then, for n > n^ 2 \ the trajectories scape 
from the influence of the period-one fixed points and I 2 
starts to increase again. 

In Fig. 3(a) we also explore the case Io ~> K, see open 
symbols. During the first few iteration steps, since K is 
small as compared to Iq , I 2 does not increase significantly 
as a function of n; so, I 2 remains approximately equal to 
Iq up to a crossover iteration number n^K For n > 
, I 2 follows the same panorama when increasing n as 
it does in the case I <C K: it grows up to , then 
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FIG. 3: (a) J 2 as a function of n in the slow diffusion regime 
(K -C 1). Full symbols (open symbols) correspond to J <C K 
(Jo 3> ^). Each curve is the average over 1000 trajectories 
having initial random phases in the interval < 9a < 2n. The 
dashed lines proportional to n and n 2 are plotted to guide the 
eye. (b) [(c)] I 2 /n 2 [I 2 /n] as a function of K for n < 
[n > ni 2) ]. The dashed line equal to 0.17K 2 [0.17K 5/2 ] is the for n > n< 2) . 
best power law-fit to the data. 



FIG. 4: (a) [(b)] [J 2 at ] as a function of K. The dashed 
line equal to I.IK^ 1 ^ 2 [0.31JC] is the best power-law fit to the 
data, (c) Scaled curves I 2 K~ X as a function of nK 1 ^ 2 in the 
slow diffusion regime (K <C I). Full symbols (open symbols) 
correspond to Jo <C K (Jo ^ K). Same data as in Fig. 3. 



for n < nil' while I 2 oc n 



it becomes approximately equal to J 2 at up to n^ 2 \ and 
finally it grows again. 

Then, based on Fig. 3(a), we postulate the following 
scaling relations: 

I 2 (n,K) oc n a K p (5) 
for n < n£f) and n > with 

n£\K)otK^ (6) 

and 

n^\K) oc ; (7) 

in addition 

ll(K)*K s . (8) 

Also, from Fig. 3(a), we concluded that — const. « 
215. Below, we present a detailed analysis that allows us 
to obtain the scaling exponents a, /?, 71,2, and S. 

By performing power-law fittings to the growth 
regimes of J 2 , we determined that a « 2 for n < n£9 
and a ss 1 when n > n[ 2 \ See dashed lines in Fig. 3(a). 
Once we know the exponents a we can extract the ex- 
ponents /3. To this end, in Figs. 3(b) and 3(c) we plot 



I 2 /n 2 for n < n£P and I 2 /n for n > n^ 2 \ respectively, as 
a function of K. The dashed lines, equal to 0.17K 2 and 
0.17K 5 / 2 , which are the best power-law fits to the data, 
prove that (5 w 2 for n < and /3 ~ 5/2 when n > n^ 2 ' . 
In fact, the dependence J 2 oc K 5 / 2 for n > n^ 2 ' is not 
surprising since theoretical results for the saw-tooth map 
[23] as well as numerical computations on the stadium 
map [9] (both discontinuous maps) show that I 2 oc K b / 2 
when K < K c for large n. More precisely, for K < K c 
the dynamics of map (2) is diffusive (after the transient 
time n£ 2 ') with diffusion rate D = D K 5 / 2 VT [8], where 
D = lim„^oo (l 2 (n)) jn and the average (• • • ) is per- 
formed over an ensemble of trajectories with the same 
initial action I and random initial phases 6*o- Here, 
D ~ 0.4 is the constant proper of the choice of f(6) 
we made [8]. 

Then, in Figs. 4(a) and 4(b) we show and J 2 at 
as a function of K, respectively. We computed as 
the intersection of a power-law fitting curve I 2 oc n 2 for 
n < with the constant curve I 2 — I 2 ,. The dashed 

cr sat 

lines in Figs. 4(a) and 4(b) equal to LIK' 1 / 2 and 0.31K, 
respectively, lead to 71 ~ —1/2 and S ~ 1. As a conse- 
quence of the scalings above, in Fig. 4(c) we present the 
scaled curves I 2 K~ X as a function of nK 1 / 2 showing the 
collapse of 7 2 t and n^K 

We want to stress that the scaling of I 2 for the dis- 
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FIG. 5: I 2 /! 2 ^ as a function of nK 3/2 in the slow diffusion 
regime (K < I). From left to right: K = 1CT 4 , 2 x 1CT 4 , 



4 x 10" 
Inset: 



10" 



2 x 10~ 3 , 4 x 10-\ and 10" 



,(2) 



1.9-ff 3 ' 2 is the best power law-fit to the data. 



o = if/100. 

as a function of K. The dashed line equal to 



continuous map of Eq. (2) in the slow diffusion regime 
obeys the same scaling laws than CSM in the regime of 
weak nonlinearity, see [13], except for the appearance of 
the second crossover iteration number rv£> . To study the 
dependence of on K, in Fig. 5 (Inset) we plot n^> vs 
K. The power- law fitting of the data leads to 72 ~ —3/2 
and a proportionality constant ss 1.9. Then, in the main 
panel of Fig. 5 we show that the curves I 2 /I^ t are prop- 
erly scaled, for n > , when plotting them as a function 
of nK z / 2 . The behavior I 2 oc n for n > nffl should be 
expected in map (2) since here, in contrast to CSM with 
K < K Cl the movement is not bounded by KAM tori 
and particles can diffuse along the phase space cylinder 
without limit. 



B. Quasilinear diffusion regime 

In Fig. 6(a) we present I 2 as a function of n in the 
quasilinear diffusion regime (K ^> 1) for several combi- 
nations of K and Iq. 

For Jo <C K, I 2 grows proportional to n for all n. See 
full symbols in Fig. 6(a). For I >• K, I 2 as a func- 
tion of n is almost constant and approximately equal to 
I 2 up to a crossover iteration number n„. Then, when 
n > n„, I 2 increases proportional to n. See open symbols 
in Fig. 6(a). This behavior for I 2 is completely equiva- 
lent to that for CSM in the strong nonlinearity regime 
[13]. That is, the scaling given in Eq. (3) is valid for 
n > n„ with a rts 1 and /? w 2. This is consistent with 
the random phase approximation [3] that predicts, for 
K > K c , diffusive motion along the I direction with a 
diffusion rate D = K 2 /2. Moreover, we observed that 
the crossover iteration number n„ scales as 

njK,I ) <xKi*q* . (9) 
To get the exponents 73^4 in the scaling relation above 




FIG. 6: (a) I 2 as a function of n in the quasilinear diffusion 
regime (K ^> 1). Open symbols (full symbols) correspond 
to J S> K (To >C K). Each curve is the average over 1000 
trajectories having initial random phases in the interval < 
60 < 2tt. The dashed lines, plotted to guide the eye, are 
proportional to n. (b) [(c)] n cr as a function of Jo [K] for 
K = 10 3 [I = 10 6 ]. The dashed line equal to 37q [S.IK ~ 2 ] is 
the best power law-fit to the data. 
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FIG. 7: Scaled curves I 2 I 2 as a function of nK 2 I 2 for 
Io S> K. The dashed line shows that I 2 oc n for n > n CI . 



in Figs. 6(b) and 6(c) we plot: (i) n„ as a function of I 
for fixed K and (ii) function of K for fixed Iq, 

respectively. Using power-law fittings, see Figs. 6(b) and 
6(c), we found that n cr cx K~ 2 Iq with a proportionality 
constant w 3. Thus, we concluded that 73 w —2 and 
74 w 2. Finally, from scaling (9), in Fig. 7 we show that 



all curves I 2 I 2 as a function of nK 2 I 2 colapse into a 
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K J2> 1 
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< n < ni 2) 






I 2 


oc ni^/ 2 


n > n< 2) 


n > ni 2) 






I 2 


oc nif 2 






n > n cr 


n > n cr 



TABLE I: Behavior of J 2 in the slow diffusion (A < 1) 

and the quasilinear diffusion (K ^> 1) regimes. We have 

found that / 2 at w 0.31A, n^ )} w 215, n^' w 1.1A' -1 '' 2 , 
n( 2 > w 1.9A- 3 / 2 , and n„ w 3I$K-' 2 . 



single one. 

III. CONCLUSIONS 

We have studied the scaling properties of the action 
variable I for the discontinuous map of Eq. (2). We focus 
on the slow diffusion (K < K c ) and the quasilinear diffu- 



sion (K > K c ) regimes, being K c = 1/T. Wc found that 
the scaling of I 2 for map (2) when K <C K c and K ^> K c 
obey the same scaling laws than CSM in the regimes of 
weak and strong nonlinearity [13], respectively. Except 
for that in the slow diffusion regime, due to the absence of 
KAM tori to bound the motion, I 2 oc n for large enough 
n. Also, we conclude that the scaling I 2 oc n a K^ applies 
to discontinuous maps with 

(i) a w 2 and /? w 2 for K <§C K c and small n; 

(ii) awl and j3 s=s 5/2 for AT <C and large n; and 

(iii) awl and (5 w 2 for A' >• AT C and large n. 
Our results are summarized in Table I. 
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